REPORT 1153 


ON THE APPLICATION OP TRANSONIC SIMILARITY RULES 
TO WINGS OF FINITE SPAN 1 

By John R..Spreitbr 


SUMMARY 

The transonic aerodynamic characteristics of wings of finite 
span are discussed from the point of view of a unified small 
disturbance theory for subsonic, transonic, and supersonic 
flom about thin wings. Critical examination is made of the 
merits of the various statements of the equations for transonic 
flow that have been proposed in the recent literature. It is 
found that one of the less widely used of these possesses con- 
siderable advantages, not only from the point of view of a priori 
theoretical considerations but also of actual comparison of 
theoretical and experimental results. The similarity rules and 
known solutions of transonic-flow theory are reviewed, and 
the asymptotic behavior of the lift, drag, and pitching-moment 
characteristics of wings of large and small aspect ratio is 
discussed. It is shown that certain methods of data presentation 
are superior for the effective display of these characteristics. 

INTRODUCTION 

The small perturbation potential theory of transonic flow 
proposed apparently independently by Oswatitsch and 
Wieghardt, Busemann and Guderley, von Kdrm&n (refs. 1 
through 6), and others is now supplying a fund of information 
regarding transonic flow about aerodynamic shapes. Solu- 
tions have been given for two-dimensional flow around 
airfoils at both subsonic and supersonic speeds in papers 
by Guderley and Yoshihara, Yincenti and Wagoner, Cole, 
Trilling, Oswatitsch, Gullstrand (refs. 7 through 16), and 
others. In the application of these results to specific ex- 
amples, two items of theoretical interest have been noted 
(see, in particular, refs. 8, 17, and 18): (a) The theoretical 
results appear to be applicable at Mach numbers far removed 
from 1 even though, in most cases, the results have been 
obtained from equations valid only in the immediate neigh- 
borhood of sonic speed, (b) In the application of theoretical 
results to specific examples at Mach numbers other than 1, 
it has been noted that certain ambiguities exist in the 
theoretical determination of aerodynamic quantities. It is 
one of the purposes of this report to investigate these two 
points. This is accomplished by examining transonic flow 
from the point of view of equations that are valid throughout 
the Mach number range rather than only in the neighborhood 
of sonic speed. Such ah approach emphasizes the relation 
between the roles of linear theory and of nonlinear theory 
in the transonic range. 

The similarity rules provided by the theory (refs. 5, 6, 
and 19 through 22) have also proved to be useful in the cor- 


relation and interpretation of experimental data. It is 
with the latter aspect of the transonic-flow problem that 
the present paper is primarily concerned. In this paper, 
the similarity rules and their application to the specific 
problem of concise presentation of lift, drag, and pitching- 
moment characteristics of wings axe given in detail. The 
known solutions of two-dimensional transonic flow are re- 
viewed and the asymptotic behavior' of the aerodynamic 
characteristics of wings of large and sinall aspect ratios is 
examined. It is shown that certain methods of data 
presentation are advantageous for displaying these charac- 
teristics. 
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1 Supersedes NAOA TN 2726 entitled “On the Application of Transonlo Similarity Rules,” by John R. Sprelter, 1952. 
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section lift coefficient 

[UMt/c)Y* (da) 
drag function 

drag function for symmetrical nonlifting wings 
drag due to lift function 

section drag function for symmetrical nonlifting 
airfoils 

section drag due to lift function 
coefficient of nonlinear term of differential 
equation for <p. (See eqs. (7), (18), (19), (20), 
(23), and (35).) 
lift function 
section lift function 
pitching-moment function 
section pitching-moment function 
local Mach number 
fre e-stream Mach number 
pressure function 

stretching factors defined in equation (B8) 
maximum thickness of wing 
free-stream velocity 

velocity components parallel and perpendicular, 
respectively, to the flow direction -ahead of a 
shock 

Cartesian coordinates where x extends in the 
direction of the free-stream velocity 
distance from wing leading edge to center of 
pressure 

ordinates of wing profiles in fractions of chord 
angle of attack 

W) 

ratio of specific heats, for air 7=1.4 
arbitrary constant 
CM»-1) 

[UoHtlc)] 2 ' 3 

(Mp 1 — 1) 

[W/e)]* rt 

ordinate-amplitude parameter 
velocity potential 
perturbation velocity potential 

SUBSCRIPTS 

values given by linear theory 
conditions at the wing surface 
conditions immediately upstream from shock 
conditions immediately downstream from shock 

FUNDAMENTAL CONCEPTS 


BASIC EQUATIONS 

The quasi-linear partial differential equation satisfied by 
the velocity potential $ of steady isentropic flow of a perfect 
inviscid gas can be expressed in the form (ref. 23, p. 25) 

2 $*$,$,,— 2$ I $ I $«=0 ’ ( 1 ) 

where the subscript notation is used to indicate differentia- 
tion andn is the local speed of sound given by the relation 


‘ a s =a 0 2 -^ Uo 1 ) 

In this latter equation U 0 and a 0 are, respectively, the ve- 
locity and speed of sound in the free stream and y is the ratio 
of specific heats (for air 7=1.4). 

Introducing the perturbation velocity potential <p, where 

<P= — U<p- S r$ (3) 

it is possible to express equation (1) in terms of the deriva- 
tives of <p as follows: 


(1 — M 0 *)<p ia + 

g[(7+ l)U oVx - |-2±i 

j[(7-l)^+^(^+^+2±i^] + 

g[(7-l)17 oPl +^(^ s +^+^ 1 V , a ]+ 


r (4) 


2^1 <Pn(Uo+<Px) + 

(Lq 

2^Hv>t(Uo-h<Px)+ 

C VQ 

n < Pyt 


If it is assumed that all perturbation velocities and perturba- 
tion velocity gradients (represented by first and second 
derivatives, respectively, of <p) are small and that only the 
first-order terms in small quantities need be retained, equa- 
tion (4) simplifies to the well-known Prandtl-Glauert equa- 
tion of linear theory 

(1 — Mo*)<Pxx+Vyv + <Pit=0 ( 6 ) 

where the free-stream velocity is directed along the positive 
x axis as shown in figure 1 and where M 0 is the Mach number 
of the free stream It is well known that equation (5) leads 


k 



x 

Figure 1. — View of wing and coordinate system. 
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to useful results in the study of subsonic and supersonic 
flows about thin wings and slender bodies but that it is 
incapable, in general, ot treating transonic flows. The 
failure of linear theory in the transonic range is evidenced 
by the calculated value of <p x growing to such magnitude 
that it can no longer be regarded as a small quantity when 
compared with U 0 . 

Second-order theory for thin wings would involve solution 
of the equation 


(1 ~M^)tp a -\-<p n -\-ip„= 


M c 


t 


7+1 

U 0 


-VxVxx- 


7—1 

■ Uo 


‘Pzi.'Pvv + ‘Pit) + JJ + 


( 6 ) 


0 ?1 — w 2 ) 2 + Oh — <Pz} 2 = On — <p *) 2 


[^-^+17^,+^] 

w) 


( 10 ) 


■ This result can be put into a more desirable form by making 
use of the following relation derived directly from equation 
( 9 ): 


77 3 n *i 2 Mp* 1 J-J 2 

Uo a 7+1 M 2 U ° 


( 11 ) 


Substitution of this expression into equation (10) and rear- 
rangement of the terms yields the desired result 


Actually, we are interested in retaining higher-order terms 
only to the extent that is necessary to allow the study of 
transonic flow. Examination of the known characteristics 
of transonic flow fields indicates that the first term on the 
right can often become of importance and should be retained. 
The remainder of the terms on the right can never become 
large for transonic flows about thin wings at small angles of 
attack and can be safely disregarded. The simplified equa- 
tion is 

(1 -1 + 2 ) 7-TT + + Pit -Iff <Px<Pxx := k t Px ( Pxi (7) 


It follows from the basic assumptions of small-disturbance 
transonic theory that equation (7) is valid everywhere in the 
flow field, both fore and aft of any shock waves, but cannot 
provide any information on the discontinuities in <p that 
represent the shock waves. This information must be ob- 
tained through consideration of the classical equations for 
the_velocities on either side of an oblique shock. Thus- 
li-Ui represents_the velocity immediately in front of the 
shock wave and Ut and V 2 represent the velocity components 
parallel and perpendicular, respectively, to U, that occur 
immediately behind the shock, the equation for the shock 
polar (ref. 23, p. 108) provides that 


v 2 2 =(Ui— u 2 y — ^ — 

7+1 


UiU 2 —a* 2 

W-UiUi+a* 3 


( 8 ) 


where a* represents the velocity of sound at a point where 
the local Mach number is unity and is commonly designated 
the critical speed of sound. It is related to the free-stream 
velocity and speed of sound by the following equation de- 
rived from equation (2): 


7+1 

2 


a* 2 =a 0 


7—1 


Uo 


(9) 


Except for the important case of the bow wave in supersonic 
flow, Z7i is not, in general, alined with the direction of the x 
axis, but is inclined a small angle. With the resolution into 
components parallel to the axes of the coordinate system and 
the retention of only the leading terms, in a manner similar 
to that used in the derivation of equation (7), equation (8) 
provides the following relations between the velocity com- 
ponents (potential gradients) immediately fore and aft of 
the shock: 


(1 —Mo 2 ) (<p x j — Ph) 2 + (<Pv l ~ <ps^ 2j r(‘Px l — 

^r+l (^+* 3 ) (f 5 +!s) fe 

( 12 ) 

This equation corresponds to the shock-polar curve for 
weak shock waves inclined at any angle between that of 
normal shock waves and that of the Mach lines. 

In addition to satisfying equations (7) and (12), the per- 
turbation potential must provide flows compatible with the 
following physical requirements: (a) The flow must be uni- 
form far ahead of the wind, and (b) the flow must be tan- 
gential to the wing surface. Therefore, the following boun- 
dary conditions are to be specified for the perturbation 
potential: 

at x=— co 

{ < Px)o = {‘Pt)<i—( , Pz)o —0 (13) 

at the wing surface W 



where bZ/bx refers to the local slope in the x direction of the 
wing surface. A systematic application of the perturbation 
analysis indicates that the boundaiy conditions specified on 
the wing surface should be simplified by approximating the 
fraction on the left by y>,/Z7 0 Furthermore, it is consistent 
with the assumption of small disturbances to satisfy this 
boundary condition on the two sides of the xy plane rather 
than on the wing surface. Equation (14) is therefore re- 
placed by 



where the shape of the wing profile is given by 

Z/c—tJ (x/c, y/b) (16) 

where f(x/c, y/b) represents the ordinate-distribution function 
and r is an ordinate-amplitude parameter. Note that, in 
general, a variation of r represents a simultaneous change of 
the thickness ratio, camber, and angle of attack. In the 
special case of a nonlifting wing having symmetrical sections, 
t is proportional to the thickness ratio ; for lifting flat-plate 
wings of vanishing thickness, r is proportional to the angle 
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of attack. In order to obtain unique and physically im- 
portant solutions, it is necessary to. assume the Eutta condi- 
tion (that the flow leaves all subsonic trailing edges smoothly) . 

Upon solving the above boundary-value problem for the 
potential, one may determine the pressure coefficient by 
means of the formula 

o„=— w.** (17) 

It should be noted that the results obtained by using the 
foregoing approximate equations might be expected to tend 
toward those of linear theory as the free-stream Mach number 
M 0 departs far from unity. This follows from the fact that 
the product <p x Vn becomes small relative to the linear terms 
under this condition. Solutions of the equations for tran- 
sonic flow found to date have all possessed this property. 

COMPARISON WITH OTHER STATEMENTS OF THE TRANSONIC-FLOW 

EQUATIONS 

As a result of minor variations in the perturbation analy- 
sis, recent papers have used at least four different relations 
for k, the coefficient of the nonlinear term in the simplified 
equation for the perturbation velocity potential. As indi- 
cated in the preceding paragraphs, straightforward develop- 
ment of the theory leads to the relation 

k=M 0 *?±± (18) 

This is sometimes simplified (e. g., refs. 22 and 24) to 



(19) 


by argu i ng that M 0 can be set equal to unity in this term 
without much loss in accuracy, since the right-hand side of 
equation (7) is merely an approximation to allow the treat- 
ment of transonic flows and rapidly diminishes in magnitude 
as M„ departs from unity. In some treatments (e. g., refs. 
16 and 21), equation (1) is divided by a 3 , and the quotient 
1/a 2 in each term is expanded in a binominal series. When 
this is done, the coefficient k of the term involving ipxVn is 


k=M 0 * 


j ~2+(r-l)Af/ j 


( 20 ) 


Still another expression for k is used by Oswatitsch in 
references 13 and 14. Two derivations are given, one based 
on mass-flow considerations (ref. 14) and the other (ref. 13) 
based on simplifying equation (1), under the assumption 
of nearly parallel flow, to 

(1— M 3 )<p xx +<p vs -)-<p u =0 (21) 

and expanding the variable coefficient 1 — M 3 in the series 


1-M*=1-M 0 S 


1 —M 3 dg 

a*— U, da; 


( 22 ) 


where M is the local Mach number and a* is the critical 
speed of sound as defined in equation (9). Comparison of 
equations (21) and (22) with equation (7) shows the coeffi- 
cient k in this approximation is 


k 


1-Mo 3 

a*— U 0 


(23) 


It should be noted that the four alternative relations for k are 
identical for M 0 =l and all but that given by equation (19) 
are zero for M 0 = 0. 

A similar situation arises in the derivation of the simplified 
equation for the shock polar. Here again the precise form 
of the expression for the coefficient k of equation (12) depends 
on the details of the perturbation analysis. The most 
important point from a practical point of view is that the 
same expression for k is used in both the equation for the 
potential and that for the shock polar, namely equations (7) 
and (12). While this point has not always been explicitly 
stated, it is actually a necessary condition for the existence 
of the transonic similarity rules. 

Although each of the above alternative forms of the 
perturbation equations has been used at least once in the 
recent theoretical investigatipns of transonic flow, the most 
widely used set of equations are those derived under the 
more restrictive assumption that all velocities are small 
perturbations around the critical speed of sound a* rather 
than around the free-stream velocity U 0 . In the latter 
scheme, the perturbation potential is defined by (see, e. g., 
ref. 6 or 20) 


ip'=— a*r+4> 


(24) 


and the resulting differential equation for <p' is 


'Y *f* 1 

¥>'*»+ -* <p'x<p't 

Uf 


(25) 


The approximate relation for the shock polar is 





The corresponding boundary conditions are specified as 
follows: 


at x = — <» 

(v’lo=Uo-a* « (1 -Mo 3 ), (<p' v )o=(<p'.)o= 0 

at the wing surface 

< 27) 


where the shape of the wing profile is still given by equation 
(16). The equation for the pressure coefficient is approxi- 
mated similarly, thus, 

(28) 

w 


This statement of the equations for transonic flow is clearly 
identical to those given previously, herein, when the free- 
stream Mach number is unity. Although the derivation of 
the a* equations requires that the free-stream Mach number 
be very close to unity, these equations have been used with 
good success by a number of authors to calculate the aero- 
dynamic forces on airfoils at Mach numbers considerably 
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removed from unity. In so doing, it has been suggested that 
it might he preferable to use more accurate relations for the 
pressure coefficient or the boundary conditions; for instance, 
it has been suggested that a* be replaced with U 0 in the 
equation for O p . This matter has been discussed at length 
in references 8, 17, and 18. Since no restriction requiring 
the Mach number to be near unity is made in the U 0 analysis, 
it is informative to examine the relation between the results 
of the a* and the U 0 analyses. —This is done ip Appendix A. 
It is found that the a* analysis, if performed in a completely 
consistent manner using equations (24) through (28), yields 
values for O p that are identical to those given by the more 
general U 0 analysis using k={y+l)/U 0 . This somewhat 
paradoxical result is achieved through the action of a number 
of compensating, effects and only applies to the pressures and' 
the forces and moments derivable therefrom. It should be 
noted, in particular, that the values of the local velocities 
and Mach numbers provided by the a* analysis for flows 
having free-stream Mach numbers other than unity are in 
error. Throughout the remainder of this report, the discus- 
sion will be based on the U 0 analysis. 

A significant case where the alternative relations for k lead 
to different results is the prediction of the critical pressure 
coefficient 0 Per , defined as the value of the pressure coefficient 
O p at a point where the local Mach number is unity. It is 
important that a reasonably good approximation be main- 
tained for the variation of 0*. with M„ because shock waves 
make their first appearance, and the airfoil first experiences 
a pressure drag when C p becomes more negative than G Pcr 
somewhere on the airfoil surface. In the present approxima- 
tion, C Pcr corresponds to that value of O p , and, hence, or 
<p x , at which equation (7) changes locally from elliptic to 
hyperbolic type. This condition is recognized by the vanish- 
ing of the coefficient of thus 

1 — M 0 2 — k (<p x ) „= 0 
or, in view of equation (17) 

0 PcT = -JJ (<Px)cr— feJJ (I Afo 2 ) (29) 

The exact relation for isentropic flow is (e. g., ref. 25, p. 28) 



The variation of 0 PcT with M 0 has been computed by use of ' 
the exact relation and each of the four approximate relations. 
The results are presented graphically in figure 2. It may be 
seen that a reasonably good approximation for 0 PcT is ob- 
tained over a wide Mach number range when k is taken as- 
given in either equation (18) or (23), and that a somewhat 
greater error is incurred when equation (20) is used. On the 
other hand, equation (19) leads to a very poor approximation 
for 0 Pcr . 

Similar comparisons can be made for local Mach numbers 
M other than unity by noting that the coefficient 1 —M*—k<p x 



of ipxx in equation (7) corresponds, in the present approxima- 
tion, to 1 — M 3 , thus 

l-M 3 =l~M 0 2 -^=l-M 0 s +^a P (31) 


The corresponding exact relation for isentropic flow is 




r -—I > 

i+V ¥ » 

r ~ 
Y-l 

-1+ 



l+^=~W 



l 1 J 

_ 


(32) 


The results so obtained are generally similar to those in- 
dicated in figure 2, although the relative accuracy of the bet- 
ter approximations changes somewhat with the situation. 
All the approximations are exact, of course, when C P — 0; on 
the other hand, none of the approximations are exact, except 
for isolated cases, when G p is different from zero, even though 
all the approximations agree among themselves when the 
free-stream Mach number is unity. In order to provide some 
information regarding the errors that are likely to be in- 
curred when G p is not very small, figure 3 has been prepared 
illustrating the variation of local Mach number with pressure 
coefficient for a free-stream Mach number of unity. 

A second case where the exact and approximate relations 
may be compared is furnished by considering the velocity 
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Figure 3. — Variation of local Mach number with pressure coefficient 

for 


jump through a shock wave. If the flow ahead of the shock 
wave is uniform and parallel to the x axis, the results may be 
conveniently represented by the shock-polar diagram in 
which Vi/ a* is plotted as a function of U t /a*. The exact 
relation is furnished by equation (8). The corresponding 
approximate relations are determined from equation (12) by 
setting <p Zv <p H , <p S2 , and <p Xl equal to zero, whereby U 0 = U u 
M 0 =Mi, and 

<=[-( 1 9,* (33) 

Once the variation of <p Zi with <p^ is determined for a given 
M 0 , the corresponding variation of Vi with Z7 2 may be 
readily determined since, for this case, 

U 2 = Vi=<p H (34) 

The variation of with for M a = 1.2 has been computed 

using both the exact and approximate relations, and the 
results are presented graphically in conventional shock- 
polar form in figure 4. This figure contains, in addition to 
curves for the four expressions for k discussed previously, a 
curve computed using 

(35) 
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Figure 4. — Exact and approximate shock polars. 


This expression for k arises in a series expansion of the right- 
hand member of equation (10). (For sake of completeness, 
the corresponding variation of 0 Vcr with M 0 is also included 
in figure 2.) Just as with the comparison of the critical 
pressure coefficient, the expressions for k given by equations 
(19) and (35) lead to poor approximations, whereas any other 
of the expressions lead to reasonably good approximations. 
A notable point is that the expression for k given by equation 
(18), that is k=M 0 2 (y-{- 1)IU 0) leads to the exact relation for 
velocity jump through a normal shock wave. 

In order to facilitate comparison with previous results and 
to achieve an economy of notation, the piesent analysis is 
carried as far as possible without specifying a particular rela- 
tion for k. That is, the equations of the analysis and reduced 
parameters with which the results are expressed are written 
containing k which may be equated to any of the five stated 
relations, or in fact, to .any coefficient that does not depend 
on x, y, z or <p. The actual values of the pressure coefficient 
and Mach number for an airfoil of specific thickness ratio, 
however, depend on which relation is selected for k. Whom- 
ever such values are given, they will be those obtained by use 
of the expression for k given in equation (18), that is, 

(18) 

The principal reason for this choice is that it appeals to 
provide a set of equations, or a mathematical model, which 
approximates certain essential features of transonic flow with 
superior accuracy. 

RELATION BETWEEN TRANSONIC THEORY AND LINEAR THEORY 

It is important to recognize that wing theory based on 
equation (7) is valid for all Mach numbers below the hypor- 
sonic range. At subsonic and supersonic speeds, equation (7) 
is of the same order of accuracy as the Prandtl-Glauort 
equation of liner theory (eq. (5)) although more difficult to 
solve. At Jlf„=l, equation (7) is identical with equation 
(25), now widely used in the study of transonic-flow problems. 

On the other hand, there is no a priori method for de- 
termining whether or not a solution of the equations of 
linear theory will be valid in the transonic range. One can 
only decide by solving the problem under the assumptions of 
linear theory and then inspecting the magnitudes of the 
terms, particularly of <p z , to see whether or not they can bo 
regarded as small quantities. _ If the terms are sufficiently 
small, the linear-theory solution is presumed valid oven 
though the Mach number may be near unity. Linearizod- 
theory solutions have been obtained for a great number of 
practical wing problems and their behavior in the transonic 
range is now well known. To review briefly: For unswept 
wings of infinite span, linear theory indicates that the 
magnitude of <p x on the surface of a given airfoil is propor- 
tional to 1/VjT— M^l ; consequently, <p x approaches infinity 
as M„ approaches unity and the theory is clearly inapplicable. 
For wings of finite span, however, the perturbation velocities 
may be large or small at sonic velocity, depending on the 
particular problem as discussed in detail in reference 26. 

I Specifically, for three-dimensional lifting surfaces of zero 
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thickness, the velocities remain finite everywhere except at 
the leading edges, their magnitudes generally increasing with 
increasing aspect ratio and angle of attack. For wings of 
nonzero thickness, however, <fx generally becomes large 
logarithmically as 1— My approaches zero;, consequently, 
linear theory is inapplicable within some Mach number 
range surrounding unity. 

Summarizing, linear theory is applicable to lifting surfaces 
of small or moderate aspect ratio at all transonic speeds, 
but fails for wings of finite thickness within a range of 
Mach numbers surrounding unity. The range of inapplica- 
bility diminishes to zero as the aspect ratio, thickness ratio, 
and angle of attack of the wing tend to zero. 

In treating transonic flows for which linear theory is 
applicable, it is often advantageous to consider the special 
case of sonic flow (M 0 =l) separately. Equation (5) for 
the perturbation potential then reduces to a particularly 
simple form 

<Pn+Vu=0 

Solutions of this equation, in conjunction with the boundary 
conditions given by equations (13) and (15), are identical to 
those of linear theory found by solving equation (5) and 
subsequently setting M 0 = 1, but can be obtained with much 
less effort. Since, in addition, the results of this simple 
theory, now generally known as slender-wing theory, are 
also applicable to low-aspect-ratio lifting surfaces throughout 
the entire Mach number range, a considerable number of 
solutions of slender-wing theory have been presented in the 
last few years (e. g., refs. 27, 28, and 29). These results 
are, of course, applicable to flows at M„=l to exactly the 
same extent as the results of linear theory. 

SIMILARITY RULES 

In reference 6, von KArm&n derived similarity rides for 
the pressure distribution, lift, drag, and pitching moment of 
airfoils in transonic flow using equations (24) through (28). 
The same equations were used in reference 20 to determine 
the transonic similarity rules for wings of finite span. The 
corresponding similarity rules of linearized subsonic and 
supersonic whig theory were also derived ' and compared 
with the transonic similarity rules in the latter reference. 
It was shown that the similarity rules of linear theory con- 
tain an arbitrary parameter and can be expressed in many 
forms, one of which coincides with the similarity rules of 
transonic flow. 

A derivation of the transonic similarity rules, based on 
the U„ equations with unspecified k, is provided in Appendix 
B. This derivation possesses the advantage of being based 
on a single statement of the problem of wing theory that 
is uniformly valid at subsonic, transonic, and supersonic 
speeds. It follows from the preceding discussion that the 
results so found are identical to those of references 6 and 
20 if k is equated to (y+l )/U 0 . The similarity rules for 
O v , C t , O m , and C D are given in Appendix B as follows: 


0p= (ujc) m ‘ T hujcTy' 3 ’ ^ l ~ M * A > 7 ' |] ( 3 6 ) 

° L= (Uok) in * [(ZT^r) 1 ' 3 ’ (3 7 ) 

~(U 0 k) U3 M [(ZJ^r) 1 ’ 3 ' ^ ~ M °~ A ~\ (3 8) 

° D = '(U 0 k) m B l(U 0 kry° 13 ’ A J (3 9) 

where the geometry of related wings is given by equation (16) : 

(Z/c)=t/(z/c, y/b) 

Equations (36) through (39) are functional equations. For 
example, equation (36) is to be interpreted as stating that 
the pressure coefficient 0 P is equal to r 213 / (Z7 Jc) 113 times some 
function T of a number of specified parameters. The 
foregoing equations have been written for flows where M 0 < 1. 
If M 0 > 1, the radical yf— should be replaced with 
■jMj— I. The functions "P, JC, -M, and T> are different, 
however, for subsonic and supersonic flow. Consequently, 
subsonic flows may be related to other subsonic flows by 
the similarity rules, but not to supersonic flows, and con- 
versely. 

ALTERNATIVE FORMS OF THE SIMILARITY RULES 

It is important to recognize that the similarity parameters 
may be combined or regrouped in any manner whatsoever, 
provided the same number of independent parameters is 
always retained. For instance, in much of what follows, it 
will be found desirable to use the square of V 1 — Af „ 2 /( UJct) 1 23 
and to replace Vl— M 0 M. with a new parameter (UJct) II3 A 
obtained by dividing Vl —M 0 2 A by V 1 —M„ 2 l(U 0 kry' 3 . In 
terms of these parameters, the similarity rules are 

Cp= (i Joky 3 T [(zjir)*' 3 ’ (UoiT)1/s A; 7 l] (40) 

° L== (l 'Joky 13 * \jUokrf 3 ’ ( Z7oiT ) 1/3 A ~] ( 41 ) 

° m= (U 0 k ) 113 ' M \_(U 0 kry /3 ’ WW 3 A ~] (42) 

° D = (Uok) V3 {juJtF 3 ’ ( Uok t),/s a ~1 (4 3) 

wherein the geometry of related wings is again given by 
equation (16). 

The similarity rules thus formulated are totally equivalent 
to those given by equations (36) through (39) but possess 
three outstanding advantages: 

(a) The indeterminacy at M 0 =l resulting from two 
parameters simultaneously vanishing is removed. 
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(b) The squaring of the first parameter avoids the neces- 
sity of changing parameters as sonic speed is passed. 

(c) The use of the parameter (UJct) 113 A rather than 

A /l -Mo A aids in distinguishing the regimes in which linear 
theory is applicable in the transonic range from those in 
which nonlinear theory must be used. Thus as ( UJct) 1 13 
A approaches zero, linear theory is always applicable, pro- 
vided, in some cases, that M 0 is not precisely equal to unity. 
On the other hand, as A becomes large, linear 

theory is not applicable in the transonic range and nonlinear 
theory must be used. 

The forms of the similarity rules given in equations (36) 
through (43) have been the source of some confusion, due 
to the multiple role that r plays in determining the thickness 
ratio, camber, and angle of attack. As can be seen from 
equation (16), all three of these geometrical quantities are 
linearly proportional to r. A more explicit statement of the 
symmetrical profiles of nonzero thickness may be obtained 
by rewriting the expression for Zjc in terms of the thickness 
ration tjc and angle of attack a rather than r, thus, 

(44) 

Comparison of equations (16) and (44) indicates that f/e 
plays a similar role to r but that it is necessaiy to introduce 
a second parameter a/ (fie). In this way, a new set of equa- 
tions expressing the similarity rules is obtained which 
correspond to equations (40) through (43), although expressed 
in terms of a and t/c rather than r. They are 


C p =T(UA,a;x/c,y/b) 

(45) 


(46) 

C m =J{(pA,Z) 

(47) 

C d =V&,2,% 

(48) 

where A, and 5 are similarity parameters 

follows: 

defined as 

[(Uok)(t/c)] 313 ’ ^ KWGAO] A > « t j c 

(49) 



( 50 ) 


Now, if the isentropic flow relation for C p in terms of M and 
M 0 




r -—I > 

i+Vw 

1 *" 
Y-l 

-1 + 

At 



1+^M> 

* J 



(51) 


is differentiated with respect to M 0 , M 0 is equated to unity, 
and equation (50) is introduced, the following relation, first 
given in reference 8 by Vincenti and Wagoner, is found: 


(52) 

Since the above-mentioned derivation is based to a certain 
extent on physical reasoning and makes use of the exact 
rather than approximate relation between pressure and Mach 
number, it is of interest from the present point of veiw to 
review the equivalent result contained in the model of tran- 
sonic flow .provided by equations (7), (12), (13), (16), and 
(17). In common with the other characteristics of transonic 
theory, the result can be expressed conveniently in the form 
of a similarity rule. Thus, recall the approximate relation 
for the local Mach number given by equation (31). 


l-M s =l-M 0 M-^<7, (31) 

Now, according to the similarity rules, transonic theoiy 
does not provide information about G v or M 0 alone, but only 

about parameters such as C v and Thus, the similarity 
rule for local Mach number is the following: 


„ M*— 1 M 0 * -l (U 0 ky* 

5 vjJoW/c)) m ~WoWic)] if3 nfiF 


C p =to~\ o', (53) 


The solutions of the equations for transonic flow obtained 
for slightly supersonic flow by Vincenti and Wagoner in 
reference 8 and for slightly subsonic flow by Cole in reference 
11 indicate that the approximate relation which corresponds 
to the Mach number freeze is the following: 



(64) 


SLOPE OF PRESSURE CURVE AT Mo-1 

Liepmann and Bryson (refs. 17 and 18) have made the 
following observations which enable the determination by 
simple and intuitive considerations of the slope of the G v 
versus M 0 curve at M 0 = 1. It is a well-known fact that, at 
slightly supersonic Mach numbers, the detached bow wave 
is far away from the airfoil and nearly normal. It is also 
well known that the Mach number downstream of weak 
normal shock is as much below unity as the Mach number 
upstream is above unity. Consequently, the Mach number 
distribution on the airfoil should be independent of Mach 
number in the neighborhood of M a = 1, that is, 


It may be recognized that this relation is equivalent to the 
exact relation given in equation (50) if k is independent of 
M 0 , as in equation (19). This is not the case, however, 
when the presently preferred relation for k, namely, equation 
(18), is used. Given equation (54), differentiation of equa- 
tion (53) shows that the slope of the pressure curve in the 
reduced parameters is 



It is a simple matter to derive the corresponding slope for 
the O p versus M 0 curve at M 0 =l, provided a specific expres- 
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sion is selected for k. If k is equated to M 0 1 (-y-{-l)/U c „ the 
final result is 


dO T \ 4_ 

7+1 



(56) 


It is interesting to note that if k is independent of M 0 , the 
slope of the pressure curve is 


dC P \ 4 

dM 0 /M c - 1 7+1 


(67) 


Thus, although equation (19) leads to the exact relation for 
the rate of change of the local Mach number, the slope of 
the pressure curve is considerably less accurate than that 
provided by equation (18) . 

The range of problems to which the foregoing results 
apply is not known at present. Inasmuch as the entire 
phenomenon appears to be connected with the presence of a 
detached bow wave which stands normal to the flow, intuitive 
considerations suggest that the results are probably appli- 
cable at least to symmetrical airfoils at zero or infinitesimal 
angles of attack hut perhaps not to airfoils at larger angles 
of attack or to wings of finite span. 


APPLICATIONS 

FUNDAMENTAL HYPOTHESES AND PRINCIPLES 

Tho remainder of this report is principally concerned with 
the deduction of the qualitative, and to some extent quantita- 
tive, characteristics of thin wings in transonic flow by .means 
of simple logical considerations based primarily on the simi- 
larity rules together with the following hypotheses: 

(a) Nonlinear theory based on equation (7) is applicable 
to all problems. 

(b) Linear theory based on equation (5) is valid for all wings 
at Mach numbers either appreciably below or above unity. 

(c) Linear theory is valid at all Mach numbers, except 
possibly very near unity, for wings of small aspect ratio. 

(d) The differential pressures between the upper and lower 
surfaces of a wing having symmetrical airfoil sections are 
proportional to the angle of attack for at least a small range 
of angles about zero. 

(e) The slope of Q f versus £„ at M 0 = 1, defined by equation 
(66), is applicable at least to symmetrical airfoil sections at 
zero or infinitesimal angles of attack. 

The consequences of the foregoing statements will be con- 
sistently pursued in the following sections in the discussion 
of tho aerodynamic characteristics of airfoils and complete 
wings. Throughout, the analysis will be restricted to wings 
having symmetrical profiles. Whenever specific results are 
to be used to illustrate the statements, they will nearly always 
bo for symmetrical-wedge or double-wedge profiles and for 
wings of triangular plan form. This choice is dictated by the 
availability of theoretical results. 

The basic principle in the following analysis is to express the 
similarity rules in such forms that the lift, pitching-moment, 
and drag coefficients can be studied for limiting values of the 
parameters with no chance for ambiguity due to indetermi- 
nate forms. In this respect, the statement of the similarity 


rules provided by equations (45) through (49) will be found 
particularly useful. 

PRESSURE DRAG OF SYMMETRICAL NONLIFTING WINGS 

The similarity rule for the pressure drag coefficient of sym- 
metrical nonlifting wings having profiles given by 

( Z/e=(t/c ) g(x/c, y/b ) (58) 

is obtained from equations (44) and (48) by setting a, the 
angle of attack, to zero 

Gd°— C I> o ='D 0 (£«,+) (59) 

where 

* 0= [(U„k) (tic)} 21 *’ (W* A ( 49 ) 


Therefore, drag results for symmetrical nonlifting wings 
should be presented by plotting the variation of the general- 
ized drag coefficient C Do with f 0 and A. 

At Mach numbers sufficiently removed from unity for lin- 
ear theory to apply, G z>0 must be independent of k since k 
does not appear in either the differential equation or bound- 
ary conditions of linear theory. This implies that the param- 
eters t,, and + must be arranged in such a manner that k is 
canceled completely from the above relation for 0 Do . The 
only alternative inside the function 'Z>„ is to form the product 
V5 Numerous possibilities exist on the outside of 'D„ 
depending on whether k is canceled by using £„, +, or some 
combination thereof. Although any of these are acceptable, 
the first will be preferred. In this way the following results 
are obtained: 


( n > ( IM/-1I ) 

' D * )l (Uoky /3 \[(Uok) (t/c)] m i 


- 1/2 ' 



K Uok) (t/c)r a]=-^L= ‘Do, (V|Mo 2 -l|A) (60) 


where it should be recalled that ‘T>„ 1 is a different function for 
subsonic and supersonic flow. Equation (60) is equivalent 
to the extended Prandtl-Glauert rule. For subsonic flow, 
D’Almbert’s paradox requires that the drag be zero; there- 
fore, for all wings, 

(<?»„) ,=0, (Sj,= o (61) 

M.<1 t.<o 


For supersonic flow, wave drag exists which depends on 
T A as well as on the plan form and airfoil section. 
The general functional relation for the drag coefficient of a 
family -of affinely related wings at zero angle of attack, as 
given by linear theory, is 


{0 Do ) 1=7==^®?^), & o ) 

M.>i -y/Mo —1 f.>o 

(62) 


Wings of infinite aspect ratio. — For wings of infinite span 
(or airfoils), equation (62), representing the functional 
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relation of linear theory for the drag coefficient, reduces to 
the following: 


K)« 

Af.>I 


m 2 


X const., 


e.>o 


const. 

= VE 


(63) 


where the value of the constant depends on the shape of the 
airfoil. Numerous experimental data show variations con- 
sistent with equation (63) at Mach numbers greater than 
about that of shock attachment. At Mach numbers closer 
to unity, however, the theoretical values provided by this 
equation are unreliable. It is evident from inspection of the 
results of linear theory itself that such a failure occurs, since 
the perturbation velocities assumed to be small in the deri- 
vation of the equations are found to become infinitely large 
as the Mach number approaches unity. . It is apparent, 
therefore, that it is necessary to resort to nonlinear theory 
for the calculation of the drag of airfoils in the transonic 
speed range. 

A similarity rule for the generalized section drag coefficient - 
of symmetrical nonlifting airfoils which is valid throughout 
the Mach number range may be obtained from equation (59) 
by setting A equal to infinity. 




.(W, - <PM U V ^ 

( t/c) m Ci ° ( t/c ) 5/3 ( _ fl <J ~ ' “»(«) (64) 


At a Mach number of unity, the similarity parameter 
vanishes and the function </„(£„) is a constant 

/■ tic , ) 5/ 3 

(ca o )—do(0)=cojist., ( Cd °) == (JJ k) 113 ^ const. (65) 

indicating that the section drag coefficients of affinely related 
airfoils are proportional to the five-thirds power of their 
thickness ratios. If hypothesis ( e ) is accepted, the variation 

of c do with £ 0 at M 0 = 1 is found to be zero for completed 
airfoils 


< 66 > 

If k is taken to be M 0 2 (y+l)JU 0 , (dCJdM p ) is given by 
equation (56) and the slope of the drag curve is 

(£t) U(4Si i 

\dM 0 J m. 


rKmrX .-^ d ( D 




(67) 


The corresponding exact relation can be determined simi- 
larly by use of equation (52) for the slope of the pressure 
curve. It is 


\dWoJM, 


-rft 


.7+1 7 + 1 


(Or) 




r- 


g<£> 


2 

7+1 


K) 




( 68 ) 


Since calculations have been made of the drag in transonic 
flow of simple symmetrical sections at zero angle of attack, 
it is not necessaiy to speculate further regarding the vari- 
ation of c do with £ 0 . At present, complete theoretical infor- 
.mation exists for the drag of symmetrical double-wedge 
airfoils throughout the transonic range. Solutions for this 
section have been obtained by transforming the nonlinear 
differential equation of the small disturbance transonic 
theory into hodograph variables and taking advantage of 
simplification of the normally difficult problems relating to 
the boundary conditions by restricting attention to polyg- 
onal profiles. The results for flows having subsonic, sonic, 
and supersonic free-stream velocities have been givon, 
respectively, by Trilling (ref. 12), Guderley and Yoshiliara 
(ref. 7), and Yincenti and Wagoner (ref. 8). The linear- 
theory solution for pure supersonic flows has been given by 
Aekeret (ref. 30). All of these results are combined on a 
single graph in figure 5. It may be seen that the preceding 



Figube 5. — Theoretical drag of double-wedge airfoil. 


remarks concerning the relation of the results of linear 
theory and nonlinear theory, and the slope of the drag curve 
at a Mach number of unity are illustrated by this comparison. 

Osawatitsch has given approximate solutions for M 0 <1 
for the pressure distribution on symmetrical biconvex 
airfoils (refs. 13 and 14) and for the pressure distribution 
and drag on NACA four-digit symmetrical airfoils (ref. 14). 
This work has recently been extended to several NACA 
6-series airfoils by Gullstrand (refs. 15 and 16). Tbeir 
drag results are generally similar to those indicated for the 
■double-wedge section in figure 5. 

In problems such as we are considering herein, the final 
test is provided by comparison with experimental results. 
Unfortunately, experimental results for the transonic speod 
range are scarce, as well as difficult to obtain, and no data 
are available for direct comparison with the theoretical 
results summarized in figure 5. At present, however, com- 
plete information, both theoretical and experimental, doos 
exist for a single-wedge section followed by a straight section 
extending far downstream. In accord with some of the 
original papers on this subject, the single-wedge section is 
considered as the front half of a symmetrical double-wedgo 
airfoil having a chord c. Solutions for this section obtained 
using transonic flow theory have been given for flows having 
subsonic, sonic, and supersonic free-stream velocities, re- 
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spcctively, by Cole (ref. 11), Guderley and Yoshihara (ref. 7), 
and Vincenti and Wagoner (ref. 8). These results are shown 
in figure 6 together with the corresponding experimental 


ly+ I) '/ 3 _ 
(// c) 5/3 Cd ° 



results of Liepmann and Bryson (refs. 17 and 18). Both 
the theoretical curve and the' experimental points are plotted 
in the same manner as presented in the original papers. As 
remarked in a preceding section, the theoretical results are 
based on a set of equations equivalent to the present equa- 
tions with k equated to (t+ 1)/Z7 0 . The small vertical lines 
on the experimental data points represent the uncertainty 
of the values. The slope of the drag curve at M 0 — 1 is no 
longer zero as indicated for the complete airfoil by equation 
(06) but takes on a positive value given originally by Liep- 
mann and Bryson (ref. 17) and readily derivable from 
equation (60). 



This figure indicates that the theoretical and experimental 
results are only in general qualitative agreement. Since it 
is shown in a preceding section that theoretical considera- 
tions suggest that a superior theory results if k is equated 
to M 0 a (7+1)/17 0 rather than (y+l )/U 0 , it is of interest to 
recalculate and roplot the present results accordingly so as 
to ascertain to what degree these thoughts are borne out by 
actual experiment. The results are shown in figure 7. It 
can be seen that the theoretical and experimental results 
are in nearly perfect agreement. Comparison of figures 6 
and 7 provides striking evidence supporting the contention 
that k should be equated to M 0 2 (y-\-l)/U 0 rather than 
(y+l)/U 0 , as has been done so often in the past. 

Wings of finite aspect ratio. — The similarity rules for wings 
of finite aspect ratio are given by equations (59) and (60). 
Although no essential simplification of the rules occurs for 
wings of small aspect ratio, the range of applicability of 
linear theory increases as the aspect ratio decreases. This 
point can be illustrated by considering the results provided 
by linear theory in a specific case. A good example to 
select for this purpose is that of the drag in supersonic flow 
of a triangular wing with symmetrical double-wedge airfoils. 
(See ref. 31.) This particular choice was made for the folio w- 

321005—55 67 


(//<r) 5/3 Cd ‘ 



ing reasons: (a) Solutions are known for all supersonic Mach 
numbers; (b) the double-wedge airfoil discussed in the pre- 
ceding sections corresponds to the limiting case of the wing 
of very great aspect ratio. The drag results provided for 
this' wing by linear theory are presented in figure 8. 



Figure 8. — Supersonic drag of triangular wing, linear theory. 


The results of figure 8 are presented in a different mannor 

in figure 9 wherein C Do is plotted as a function of £, for 
various values of A as suggested by equation (59). For 
purposes of comparison, the curves for the drag of airfoils 
(1= CO ) computed by both linear and. nonlinear theory are 
also included in the graph. As noted in the preceding 
section on airfoils, comparison of the results of linear theory 
with those of nonlinear theory for wings of A= <*> shows 
that good agreement exists for larger ? 0 , but that at smaller 
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fa 

Figure 9. — Variation with. {„ of reduced drag coefficient of triangular 

wing. 

the values predicted .by linear theory become too large. 
For wings of finite A, only the results of linear theory are 
available. They also exhibit the trend of indicating infinite 
drag as £„ approaches zero; however, the range of £, in 

which the value of C Do is excessive is much less than is the 

case for A— . In general, 0 Do of wings of small aspect 
ratio diminish es with decreasing aspect ratio. 

The drag results of figure 8 are presented in still another 

form in figure 10 in which is plotted the variation of C Do 
- with A for various values of The principal merit of 
this method of plotting is that it aids in distinguishing the 
region where nonlinear theory must be used from that where 
linear theory may be useful. Thus, the two-dimensional 
nonlinear theory results appear on the right of the graph 
corresponding to large A; whereas the three-dimensional 
linear theory results appear on the left for small A. The 
filling in of the remainder of the graph requires either the 
solution of the equations of three-dimensional nonlinear 
wing theory or the use of experimental data. It should 
again be noted that the present drag considerations apply 
only to the pressure drag. Before plotting experimental 
results in the manner indicated, it is necessary to first, 
subtract the friction drag. 

LIFT 

The similarity rule for the lift coefficient C L of a family 
of wings having ordinates Z given by equation (44) 

Z/c=lfc [s(*/c, y/b)- ± f] (44) 

is given in equation (46) as 

<?*=£(*., Z, 5) ' (46) 



Figure 10. — Variation with A of reduced drag coefficient of triangular 

wings. 


where 

f M,*-l 
[UoKm m> 


A=[U 0 k(t/c)] m A, 3=£ 


(49) 


If hypothesis (d) is accepted, C L varies linearly with a for 
at least small angles of attack. It is advantageous, therefore, 
to consider the lift ratio Oja rather than C L alone, thereby 
minimiz ing the influence of 3. 

1T = « (ZjJ) 1 ' 3 ' C(fo ’ 3)=S [C7 e fc(</c)] 1/ » ^’ 5) (?0) 

where £' (£,,, A, 3) is a new function of the indicated variables 
obtained from £ (£„, A, 3) by division by 3. Therefore, 
lift results may be presented by plotting the variation with 

£„, A, and 3 of a generalized lift ratio Oja defined as follows 
(the prime on £ has been omitted for simplicity): . 


C L /a=[UMm m (C L /a)=£(U A, 3) (71) 

Equation (71) shows that CJa depends upon three 
parameters, one more than the number which can readily bo 
treated on a simple plot. Simplification can be gained, 
of course, by holding one of the parameters constant for an 
entire graph. Results so presented aro particularly interest- 
ing for £<,= 0, (M 0 = 1); A= , (A=<»); and 3=0, (a=0). 
The latter scheme is especially good since experiments 
indicate that lift curves of wings are often relatively straight 

lines at all Mach numbers. The values of Oja at 3=0 
might, therefore, be expected to be good indications of the 
actual values for other 3. The appropriate similarity rule 
may then be written 

(CJ«)=£tio, A, 0)=£ o (U A) 

a =0 


(72) 
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For cases where linear theory applies (hypotheses (b) and 
(c)), two statements can be made immediately which 
provide further information about (a) (C L /a)i must 

bo independent of k, and (b) (C L /a)i must be independent of 
3 by virtue of the superposition principle of linear theory. 
Therefore, following the procedure used in equation (60) 
gives 

(vim/- 1 1 a ), {^) r ^ (VIEW 

(73) 

where again JCi is a different function for subsonic and super- 
sonic flow. 

Wings of infinite aspect ratio. — For wings of infinite aspect 
ratio, the functional relation of linear theory for the lift ratio 
given by equation (73) reduces to 




const. 

VW^TI 


©- 


const. 


Solutions of the equations of linear theory show that the 
value of the constant is 2ir for subsonic flow and four for 
supersonic flow. Examination of these results indicates that 
they are valid at Mach numbers appreciably less than or 
greater than unity, but are invalid for Mach numbers near 1. 

A similarity rule for the section lift coefficients of a family 
of affinely related symmetrical airfoils which is valid through- 
out the Mach number range may be obtained from equation 
(71) by setting A— <=° . 

(C[/a)=-C(£o> 00 j «) =(£<>, la) (75) 

At a Mach number of unity, ( 0 is identically zero and the 
expression for the lift ratio becomes 

Equation (76), when considered together with hypothesis (d), 
indicates that at sonic speed, the lift-curve slope at zero 
angle of attack of airfoils of a single family varies inversely 
as the cube root of the thickness ratio, thus, 

(c ,/a) %o 1 = [uM/cW 3 /(0, a ^ Q)= [Uokm\ m (77) 

Note that as the thickness ratio goes to zero, the value of the 
lift-curve slope at zero angle of attack becomes infinite, just 
as is indicated by equation (74) to be the case according to 
linear theory. If, on the other hand, t/c diminishes to zero 
while a is fixed so that 3 becomes very large, it is plausible 
that the thickness ratio does not have any effect on Cj. We 
thus have the following: 

^^°0 = ((w /3Xconst - (78) 

In this case, c, is proportional to the two-thirds power of a. 

If hypothesis (e) is acceptable, the variation of cja with 
at M 0 =l is zero, since 


[k flw. (I), X<!) 

-SKtLXC)- 0 ' cm 

If k is taken to be M 0 2 (y+i)/U 0 , (dG,fdM 0 )it.mi is given by 
equation (56) and the variation of Cj/a with M 0 at M 0 = 1 is 


[dMo C'LL-rf O-fUU 

M.= l 

1 (f) f(SL (80) 

Alo«=l 

The corresponding exact relation can be determined similarly 
when equation (52) is used for the slope of the pressure curve. 

[dll "TFT GL (81) 

AT. -I 

At present, calculations have been made of the lift in sonic 
flow (Guderley and Yoshihara, ref. 32) and in slightly super- 
sonic flow (Vincenti and Wagoner, ref. 10) of symmetrical 
double-wedge airfoils inclined an infinitesimal angle of attack. 
The corresponding solutions for airfoils in slightly subsonic 
flow have not yet been found. The results of the above- 
mentioned investigations are shown in figure 11 together with 


Prondtl- Glauert- V 



Figube II. — Theoretical lift-curve slope of double-wedge airfoils. 

the corresponding values given by linear theory. As in the 
drag case discussed in the preceding section, it may be seen 
that the remarks concerning the relation between the results 
of linear theory and nonlinear theory and the slope of the 

curve of c,/a versus £„ at f o =0 are verified for this particular 
airfoil. 

Wings of vanishing aspect ratio. — Two well-known results 
of linear theory are that the lift-curve slopes of wings of 
finite aspect ratio remain finite throughout the entire Mach 
numbor range and that the lift-curve slopes of wings of 
vanishing aspect ratio are independent of Mach number. 

Therefore, equation (73) implies that ( C L /a) t must be pro- 
portional to A either for wings of vanishing A in any flow, 
or for any wing in a flow of vanishing 

(GjJ a) i = (OlI a)i=Ax const. , ifJhl a)i—(OjJ «);= A Xconst. (82) 

A-*0 f.~0 A-+0 

The value of the constant must be determined for each plan 
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form by actually solving the equations of linear theory. 
For wings having trailing edges which possess no cutouts 
extending forward of the most forward station of maximum 
span, that is, triangular, rectangular, elliptical, etc., as well 
as certain sweptback wings, the value of the constant is 
tt/2. (Refs. 26 through 29 should be consulted for further 
discussion of this point as well as for the values of the con- 
stant for -wings having cutouts in the trailing edge which 
violate the above-stated condition.) 

It is seen from equation (82) that the lift-curve slopes of 
wings in sonic flow decrease continuously in magnitude as 
the aspect ratio diminishes toward zero. Since, in addition, 
the lift-curve slope given by linear theory has its maximum 
value at M 0 — 1, it is conjectured that the lift results of linear 
theory are a good approximation to those of nonlinear 
theory not only for all wings at Mach numbers far from unity 
but also for all Mach numbers for wings of sufficiently small 
aspect ratio. 

Wings of finite aspect ratio. — At the present time no solu- 
tions of the nonlinear theory are available for wings of finite 
aspect ratio. However, from the remarks of the preceding 
paragraphs, it is apparent that a curve representing the varia- 
tion of CJa with A for constant £„ and a would have the 

following asymptotic properties: C T Ja. would increase linearly 
with A for small A and be independent of A for large A. 
In order to give a better idea of the numerical values to. be 
expected, a set of typical results of this type is shown in 
figure 12 for wings having triangular plan forms and sym- 



Figube 12. — Variation with 2i of reduced lift-curve slope of triangular 

wings. 


An interesting result of the foregoing remarks concerns the 
influence of thickness ratio on the lift-curve slope at zero 
angle of attack of wings in sonic flow. For wings of largo 
aspect ratio, the lift-curve slope is inversely proportional to 
the cube root of the thickness ratio. For wings of small 
aspect ratio, the lift-curve slope is independent of the thick- 
ness ratio. 

PITCHING MOMENT 

The remarks of the pitching-moment characteristics of 
wings follow in a manner exactly analogous to those just 
stated for the lift characteristics. The corresponding state- 
ments for the pitching-moment coefficient C„ may be obtained 
by simply replacing C L with C m and £ with M. Thus, the 
similarity rules for C m corresponding to equations (71) and 
(73) are the following, respectively: 

(CM= [U 0 km] m (O m /a)= MS', A, 5) (83) 

(84) 

where once more Jtf and M x are different functions for sub- 
sonic and supersonic flow. The only difference between tho 
discussion of C m and 0 L is that the values of the constants of 
equations (74), (77), (78), and (82) are, of course, different. 
Graphs of theoretical pitching-moment characteristics for 
airfoils and for triangular wings corresponding to tho lift 
results of figures 11 and 12 are shown in figures 13 and 14. 


■ 
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Figure 13. — Theoretical pitching-moment curve slope of double-wedge 

airfoils. 


The moment axis is taken to be through tho most forward 
point of the wing. 

Sometimes it is deshed to present pitching-moment char- 
acteristics of wings in terms of center-of-pressure position 
rather than pitching-moment coefficient. Since the centor- 
of-pressjire position can be expressed in terms of O m and G h by 


* c.p . G m 
c ~ C L 


( 86 ) 


metrical double-wedge airfoil sections. The supersonic 
results are those of Stewart, Brown (refs. 33 and 34), and 
others. The subsonic results are those calculated by De 
Young and Harper (ref. 35) using Weisslnger’s modified 
lifting-line theory. 


the resulting expression for the center-of-pressure position 
found through application of equations (71) and (83) is 

%c.p. -A? oj) 

mo, A, a) 


c 




( 86 ) 
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Figure 14. — Variation with A of reduced pitching-moment curve slope 
of triangular wings. 

Tlie corresponding relation for linear theory is 

(~t)~C V-i (JW7=l\A)=C 'P-.CVf&l A) (87) 

PRESSURE DRAG DUE TO LIFT 

The similarity rule for the pressure drag of inclined wings ' 
having symmetrical airfoils is indicated by equation (48) to 
be the following: 

K= ( jj^C D ='T>(l; 0 ,A,Z) ( 88 ) 

The portion of the drag due to lift is therefore 

A C D =0 D -0 D =^j^m^ A, 3)-m, A, 0)] 

A, a) - (89) 

Since the differential pressures between upper and lower sur- 
faces of wings having symmetrical airfoil sections are pro- 
portional to a for at least a small range of a surrounding zero 
(hypothesis (d)), it follows that AC' d is proportional to the 
square of a. It is advantageous, therefore, to consider the 
drag-rise ratio A Cd!<A rather than A C D alone, thus, 

^ < ‘ 90) 

where 'Da' is a new function of the indicated variables ob- 
tained from 'Da by dividing by 3 s . Consequently, drag-due- 
to-hft data should be presented by plotting the variation with 

f 0 , A, and 3 of a generalized drag-rise ratio AO D jo? defined as 
follows (the prime on *Da is omitted for simplicity) : 

AO^^[U 0 h(tlc)Yi\AC D lcA)=^^oX S) (91) 


The actual presentation of the results of this three-parameter 
system may be accomplished as described in the section on 
the lift of wings. ' Of particular interest is the simplification 
resulting from presenting only the values found at a= 0. 
The simplified similarity rule is then 

(A^M='Da(£,,2,0) (92) 

a-o 


For cases where linear theorv applies, the following results 
hold: 



1 






(93) 


Wings of infinite aspect ratio. — For wings of infinite aspect 
ratio, the functional relation of linear theory for the drag due 
to lift, equation (93), reduces to 


AcA const. /AcA const. 

. A~y|M 0 2 -ii’ v ) r yf^ 


(94) 


Solutions of the equations of linear theory show that the 
value of the constant is zero for subsonic flow and four for 
supersonic flow about any symmetrical airfoil. These results 
are valid at Mach numbers appreciably less than or greater 
than unity but are invalid for Mach numbers near- 1. 

A similarity rule for the drag due to lift of a family of 
affinely related symmetrical airfoils that is valid throughout 
the Mach number range may be obtained from equation (91) 
by setting A= co . 

» , S)= &(&, S) (9 5) 


At a Mach number of unity, equation (95), for the drag due 
to lift, reduces to the following: 

( (~?) = [ ujmF dA ( 0> w) (96) 

f.=0 M.-l 

Equation (96), together with hypothesis (d), indicates that, 
at sonic speed, the drag-rise ratio A c d /a? of airfoils of a single 
family varies inversely as the cube root of the thickness ratio. 
For very large values of 3, the thickness ratio cannot have 
any effect on Ac d ; therefore, A c d is proportional to the five- 
thirds power of the angle of attack. If hypothesis (e) is 

accepted, the variation of Ac d /a? with £<, at M 0 =l is zero, 
since 


-f§( < 97) 

If k is taken to be My (y+l)/U 0 , (dC p /dAI 0 )u 0 -i isTgiven by 
equation (56), and the variation of Acj/a 2 with M 0 at M 0 =l is 
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n* i mM 

\jiMo\oc J a 


dx |_ dM 0 


! L/(?) 


=^i[^i)-|(v)J<?)“-!(vL (98 > 

The corresponding exact relation can be determined similarly 
by use of equation (52) for the slope of the pressure curve. 

[dir,, G?) (99) 


a-K) 


Wings of vanishing aspect ratio. — Since, as'indicated in 
the preceding sections, the lift calculated by means of linear 
theory remains finite throughout the Mach number range, 
and di mini shes to zero as the aspect ratio approaches zero, 
it follows that the drag due to lift must behave s imil arly. 
It is consequently reasonable to presume that linear theory is 
capable of describing the drag-due-to-lift characteristics of 
wings of vanishing aspect ratio at all Mach numbers. 
Therefore, the following relations stemming from equation. 
(93) hold: 


(ix), “(tt). 

a- k) r.-o 

( ’tt ), -(^). -^ Xco ““- 


A— 0 JWo*=J 


( 100 ) 


Solutions of the equations of linear theory show that the 
value of the constant is tt/ 4 for all wings of small V|M 0 S — l|A 


whose trailing edges possess no cutouts extending forward of 
the most forward station of maximum span (i. e., triangular, 
rectangular, elliptical wings, etc.) . The quoted value of the 
constant corresponds to the development of the full “leading- 
edge force.” It is known .that this force is oftentimes not 
completely realized, due to a local separation and subsequent 
reattachment of the flow around the leading edge. If the 
leading-edge force is nonexistent, the corresponding value 
for the constant is v/2. 

Wings of finite aspect ratio. — At the present time no drag- 
due-to-lift results have been obtained from the nonlinear 
theory for wings of either finite or infinite aspect ratio. 
The foregoing remarks, however, are sufficient to determine 

that a curve representing the variation of A Od/cc* with A 
for constant £ 0 and a would increase linearly with A for 
small A (unless the degree of attainment of the leading-edge 

zv <■</ 

force also depends on A .) and become independent of A for 
large A. The resulting curve would presumably have the 

same general appearance as that shown in figure 12 for CJa. 

It may sometimes be desired to present drag-due-to-lift 
results in terms of A C D /C L 3 or AG D laO L rather than A C D jo?. 
The similarity rules for these quantities can be quickly de- 
duced from the foregoing results. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Dec. 16, 1962 



APPENDIX A 


RELATION BETWEEN U Q AND a* STATEMENTS OP THE TRANSONIC-FLOW EQUATIONS 


Equations (3), (7), (12), (13), (15), and (17) were presented 
in the text as being applicable to the study of transonic, as 
well as supersonic, flow about thin wings. These equations 
repeated below as equations (Al) through (A6), will be 
referred to as the U 0 statement of the problem since the 
perturbation velocities are taken about the free-stream veloc- 
ity U 0 . The perturbation potential <p is defined by 

<P= — U l #+$ (Al) 


at the wing surface 

(?',),-o=a*(|D (All) 

and the pressure coefficient is given by 

0 T =~ Je [<p'x—(<p'k)o\ (A12) 


The differential equation is 

(1— M 0 2 ) <p«+ <Pn+<Pz,=k<p x <p xx (A2) 

/ 

The shock relation is 


(l—M/)(<p x — <p H ) 2 +(<p tl —<Pv i y+(‘Pz 1 — <p l2 ) 2 


-*(^“)k -.)■ 

(A3) 

The boundary conditions are 
at x=— co 

(<Px)o OPv) 0 (<Cj)o 0 

(A4) 

at the wing surface 


' M— G§) 

(AS) 

The pressure coefficient is given by 

& 

1 

II 

(A 6) 


The relation between the U 0 and a* statements can be 
determined directly in the following manner. The differ- 
ential equations and shock relations for ip' and <p will be the 
same if 


T+l V 
a* bx 



-(1 -Mo 2 ) 


or if 


ip' — k 


a* 

7+1 



(1 -M*) 

7+1 


x 


(A13) 


The boundary conditions for tp' corresponding to those 
stated for <p in equations (A4) and (A5) are: 
at x—~ co 


W*)o=k 
Op v)o~k 


a* , s a* (1—M 2 ) a*(l-M„ s n 

7 + 1 7+1 ~ 7+1 

W«= o, {<p',)o= k {<?,)„= o 


(A14) 


at the' wing surface 


In the usual a* statement of the transonic-flow equations 
(eqs. (24) through (28) in the text), it is assumed that all 
velocities are only slightly different from the critical speed 
of sound a*. The perturbation potential is defined by 


(<p'z%-o—k (<Pz)z-o—kU 0 yq-y (Alo) 

Finally, the pressure coefficient is given by 


<p'=—a*x+$ (A7) 

The differential equation is 

<p'x<p'tx (A8) 

l v 

The shock relation is 

(*-/-) (A- A) 2 (A9) 

If the perturbation analysis is carried out in a completely 
consistent manner, the boundary conditions are: 
at x= — oo 

>g.=(iW- D~v w,=w t = 0 (Aio) 


C»'=-| ^ 0, (A16) 

Comparison of the above equations with those given 
previously for the completely consistent a* analysis reveals 
their identity if, and only if, k is taken as (y+l)/U 0 , that is, 

O v '= C r if (A17) 

U 0 

As far as obtaining the values of 0 P is concerned, therefore, 
the conventional a* statement of the problem may be 
regarded as a transformation of the more general U 0 state- 
ment with k being defined as in equation (A17) or (19). As 
is evident from equation (Al3), however, tbe local velocities 
found in the a* analysis are only correct when the free- 
stream Mach number is unity. 
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Tlio foregoing discussion has demonstrated that the con- 
ventional a* statement of transonic-flow problems yields 
the same values for C p as the U 0 statement, provided k in 
the latter set of equations is equated to (y-\-l)/U 0 . It is 
indicated in the text, however, that transonic results ob- 
tained by use of other expressions for k (and, in particular, 
Jc=Mf(y+l)IU 0 ) are superior in many important aspects 
to those obtained using k=(y+l)/U 0 . It is consequently 
not without interest to determine, a generalized form of the 
a* equations which correspond to the U 0 equations with 
unspecified k. The resulting equations for the perturbation 
potential and the shock relation are: 


at the wing surface 

-o= a* (|D (^22) 

and the pressure coefficient is given by 

<V=-4 (A23) 

Ur 

It can be readily verified that these equations will produce 
the same values for the pressure coefficient as equations (Al) 
through (A6), provided k and /:' are related as follows: 


<p'.=— a*x-\-$> 

<p'vi>-\~ , p'ii == k'<p'x<p'xx 


(Al8) 

(A19) 


k a v 


(A24) 


Wv-'P , v 2 Y+{<P , '-'p'h) <1 =V (A20) 

The boundary conditions are: 

atr=— oo 

(<p'v)o=(<p',)o=0 


Thus, if it is desired to obtain values for O p by use of equa- 
tions (Al8) through (A23) that are the same as those given 
by equations (Al) through (A6) with k=M 0 a (y+\)/U 0) k' 
should be equated, to the following: 1 

t , Mp (t+1) 
a* 


(A21) 


(A25) 



APPENDIX B 

DERIVATION OF SIMILARITY RULES 


The basic equations of linear theory and of nonlinear theory 
of transonic flow may be summarized as follows. The 
differential equations are: 

/ 

|0 Linear (Bl) 




I k ~ — ’ Nonlinear 


(B2) 


The approximate shock relation used in the treatment of 
transonic flow is 


( 1 — M o 2 ) ( ^ ) 2 + ( <p H — <pj) s - + ( <p, l — ) s 


The boundary conditions are: 
atx= — co 


=i .(SL+fs) («,-«,)■ 

i: 


at the wing surface 

J_ d f A 2/N 

17. W.-o b(x/c ) 1 \c b) 

where the geometry of the wing is given by 
(Zlc)=rf(xlc, y/b ) 

The pressure coefficient is given by 

p 2 by 

U 0 bx 


(B3) 


(B4) 


(B5) 


(B6) 


(B7) 


If the differential equations are now transformed into a 
system with primed quantities and the proportionality or 
stretching factors are denoted by s with appropriate sub- 
scripts such that 


x'==s x x, y'=s r y, z'=8,z, <p'=s 9 <p, TJo —&u Uo 

/(l -'I??/ Vl ~ M o\ U 0 'k'=8 t U 0 Jc or k'=^k 
the potential equations (eqs. (Bl) and (B2)) become 

(l _ M av, J2 

1VJ.0 ) >u,/jT „ 5\„/2 


(B8) 





Linear 


(B9) 


Nonlinear (BIO) 
StSp 2 bx' bx' 2 s 2 


and the shock relation becomes 




S$£ 


8 / 




_ 8y8 x 

"s*V 


Nonlinear (Bll) 


The flow in the primed system is similar to that in the 
original system if y' satisfies the same differential equations 
and boundary conditions as y. Consequently, for similarity 
to exist, it is necessary first of all that the potential equations 
and the shock equations for the two flows be the same. 
Therefore, the following relations between the stretching 
factors must be satisfied: 

Me-i £»_i 8 - ( X/X ' hhieax - (B12) 


s x 


8x8u 


Cx/x' 

W/«* 


Nonlinear 


(B13) 


where, for linear theory, A/X' is an arbitrary constant which 
can be equated to «//«* if desired. The constant is written 
as a fraction in order to maintain a certain symmetry 
throughout the analysis. 

An immediate consequence of this transformation is that 
the wing plan forms undergo an affine transformation such 
that the aspect ratios of -wings in similar flow fields are 
related, according to both linear and nonlinear theory, by 


or by 


At £? A —A / ^ Mo J 

A ~8x 8 „ A ~i (1 —Me 2 )' A 

V (1 -Mo 2 )' A'= Vl -M a * A 


(B14) 


(Bl 5) 


Since y' is proportional to y, the boimdary conditions at 
x=— co are automatically satisfied. The boimdary condi- 
tions at the wings may be given in either of two forms: 


fby'\ = s 2 /by^ = s, 
V.'-o s, V 


\bz 


UoT 


AVi) 


(B16) 


^b 2 /,_ o ° b(x/c)- 

(B17) 

whence, if the two wings have the same ordinate-distribu- 
tion functions, that is, if f(x'/c', y'/b')=f(x/c, y/b), the 
ordinate-amplitude parameters are related as follows: 


Sz;8 m $U$x 


x Ki-M*y 

X'A 1 -MS 


Linear 


(BIS) 




r (nw/ ujc\ 
~L 1 -M/ J KUo'k'J 1 


Nonlinear (Bl9) 
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or as 


\'t' Xt 


Linear 


(B20). 


A /(l —Mo 1 )' VI- -M? 

{Uo'h'rJ^iUokT)^ 


Nonlinear (B21) 


The relation between the pressure coefficients at corre- 
sponding points on the wing surface is given by 

t 2 ( dip' \ / 2 

p = U7 WA'=o _ SoS* V ~ 


Uo^xJxmO 

Linear 


(■^jr C v Linear (B22) 

NonUnear <“">> 

or more completely by 

c; vci-i ky a’-, V |i] 

p ’ MEw, i(i 7u iy a . y ~| 

Cp L(C7 0 'X*'r') 1/3 v( Mo) A ’ C' b'J 


Linear 


(B24) 


\fT^W} 

Wokr) 11 *’ 


Vl —Mo A] . Nonlinear (B25) 


The foregoing relationships may be summarized in 'the 
following statement: The similarity rule for the pressure 
coefficients on a family of wings having their geometry given 
by 

(Z/c)=t/(x/c, y/b) (B26) 

is 

Ti ( V ^%V1~MV A; f, Linear (B27) 


°'-^ y [(S^' Vi3 ^ 4 f’f] Nonlin ' or 03281 

The similarity rules for the lift, pitching-moment, and drag 
coefficients given by the linear theory are 


A) 

Linear 

(B29) 

Au vi -ms a) 

Linear 

(B30) 

<7n«=L D, ft 1 —*'* Vl -Mo 2 A) 

Linear 

(B31) 


The corresponding similarity rules given by nonlinear theory 
are 


t 2/3 

Cl= (UM' 3 * 


j 


( Uokr W 


i 


j) -yjl—Mo A Nonlinear (B32) 


—2/3 r/7 —M 2 , 1 

0 m= {u j c) m M L }Wc^T n > ^ 1 ~ Mo 2 A ] NonIineftr C 833 ) 

*2) -4] Nonlinear (B34) 

It should be noted that the foregoing equations have been 
written for subsonic flow where M„< 1. If M 0 > 1, the 
radical Vl — Af 0 2 should be replaced with ^M 0 2 — 1. 

In the linear-theory analysis, X has remained a completely 
arbitrary coefficient to be selected as best suits the particular 
problem at hand. For instance, the compressible-flow rela- 
tionships between two wings having identical pressure dis- 
tributions are found by setting X=l. If, on the other hand, 
X is set equal to Vl —M„ 3 , the thickness ratio, camber, and 
angle of attack of related wings are identical. The greatest 
simplification of the similarity rules of linear theory occurs 
when X is set equal to Vl ~—M}{r since the number of param- 
eters necessary to show the results of linear theory is thereby 
decreased by one. This degree of arbitrariness in the 
similarity rules for linear theory is in contrast to the case for 
nonlinear transonic theory in which no undetermined co- 
efficient like X is to be found. 

For the present purpose of gaining a better understanding 
of transonic flows, the most advantageous choice for X is 


(U„k) 1/3 

t V3 


(B35) 


because then the similarity rules for linear theory assume 
forms identical to those for nonlinear transonic theoiy. This 
is important since it implies that solutions of linear theory 
and of nonlinear transonic theory can be expressed' as func- 
tions of the same parameters and, hence, can both be plotted 
on a single graph in terms of one set of parameters. The two 
theories would, of course, yield two distinct curves on such a 
plot. The curve for linear theory would be accepted as valid 
for purely subsonic and purely supersonic flows but may or 
may not be valid in the transonic range, as discussed pre- 
viously. The curve for nonlinear transonic theory is valid 
not only for transonic flows, but for subsonic and supersonic 
small perturbation flows as well. As is evident from the 
derivation of the basic equations, however, the results of tho 
nonlinear'transonic theory should be considered to be of only 
equal accuracy to those of linear theory in the definitely 
subsonic and supersonic regimes, despite the fact that the 
solutions are much more difficult to obtain. 
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